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' Abstract 

The interaction of an open system S with a pre- and post-selected environ- 

c , 

, ment is studied. In general, under such circumstances S can not be described 

^3 . 

in terms of a density matrix, even when S in not post-selected. However, a 
^ ! simple description in terms of a two-state (TS) is always available. The 

I two-state of S evolves in ti me from an initially 'pure' TS to a 'mixed' TS 

and back to a final 'pure' TS. This generic process is governed by a modified 
Liouville equation, which is derived. For a sub-class of observables, which 
can still be described by an ordinary density matrix, this evolution gener- 
ates recoherence to a final pure state. In some cases post-selection can even 
suppress any decoherence. 



^e-mail: reznik@physics.ubc.ca 



1 Introduction 



The interaction of an open quantum system with an environment 0] is traditionally ana- 
lyzed while assuming a given, not necessarily known, initial state of the total closed system. 
In this case the open system can be described by a reduced density matrix which is ob- 
tained by tracing over the unknown environment's degrees of freedom. In this work we 
investigate circumstances in which the environment, and possibly also the open system, are 
bound to satisfy, not only an initial condition, but also a second final condition. In other 
words, we shall consider the interaction of a pre- and post selected environment with an 
open system. Although, under usual circumstances, such a post selection is not realized, 
it is in principle not forbidden. Quantum Mechanics is (dynamically) time symmetric, and 
it is possible to conceive situations in which the initial and final conditions are selected 
according to some 'dynamical principle' (e.g. 0). 

We shall show that when the environment is post selected, the system can not generally 
be described in terms of a reduced density matrix. At any intermediate time, between 
the pre- and post-selection, there exists no pure or mixed state, which yields the correct 
probabilities for measurements in the open system, even when the open system is not post- 
selected. We suggest that is such cases, it is preferable, both practically and conceptually, 
to describe the open system by a new object which is a generalization of the density matrix. 

It was recently suggested, that a quantum system should basically be described by 
an extension of the ordinary quantum state (or density matrix) called a "two-state" (TS), 
which is determined by two, initial and final, conditions |^ Q . In the following we apply the 
formalism developed i n Ref. |^ to this problem. The probabilities for any measurement 
in the open system are shown to be derived from a reduced TS, i.e. the TS obtained by 
tracing over the environment's degrees of freedom. When the initial and final state of the 
environment are given by pure states, this reduced TS evolves in time from an initial 'pure 
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TS' to a a 'mixed TS' (of entangled form) at intermediate times, and finally back to a pure 
TS. Therefore, the effect of post-selecting the environment is to "recohere the TS". This 
process is dynamically expressed by a modified Liouville equation. As we shall show, the 
coefficients of the new terms in the this equation are time dependent, and tuned in such a 
way that the TS finally "recoheres". 

It is well known, that interaction with an environment often causes decoherence in the 
open system. (For example see: ||, 0, ||, |, [1^). In our case of post-selection, although 
the description in terms of a (pure or mixed) density matrix is generally invalidated, one 
can still find an effective density matrix for a limited class of observables. We show that 



the post-selection causes this effective density matrix to recohere to a final pure state ITT 



In some cases, depending on the nature of the interaction, post-selection can suppress any 
decoherence. 

This article proceeds as follows. In the next section we review shortly the two-state 
formalism of quantum mechanics and elaborate on some relevant details. In Section 3. 
we apply this formalism to the case of a pre and post selected system. A simple solvable 
example is given in Section 4. The modified Liouville equation which is satisfied by the 
two-state is derived in the last section using a perturbative approximation scheme and is 
applied to some cases. In the following we set h= 1. 



2 Quantum mechanics in terms of two-states 

Two-states are particularly suitable in situations with two or more condition on a sin- 
gle quantum system. We now briefly review this formalism following Reference and 
elaborate further on some relevant issues. 

Consider a system S with a given Hamiltonian Hg. Let us assume that at ti and t2 
a complete set of measurements determine the states of S to be |'?/'m(ti)) and \ipoutit2)), 
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respectively. Now consider an ensemble of such identical systems which is defined by the 
latter two conditions. We are interested in probability distributions of observables that 
are measured in some intermediate time t2 > t > ti. The peculiarity of such a situation is 
that in general (as we shall see) these probabilities can not be derived from a single wave 
function or density matrix. It was therefore, suggested that the "state" of S at intermediate 
times should be described by a generalization of the ordinary wave function, which we call 
a 'two-state'. Generically, a TS, which we denote by g, is a non-Hermitian operator with 
the form: 

g = |...)(...| (1) 

At the left and right slots of g one inserts the information due to the conditions at the ti 
and t2 respectively. In the case of a closed system S we have: 

g{t) =U{t- ti)\^in){'4^out\U\t2 -t) = \iJin{t)){iJout{t)\ (2) 

where U (t) is the unitary evolution operator. 

More generally, two-states are elements of a Hilbert space Hn, which is defined as 
follows. Given by a Hilbert space of states Tij = we can construct the linear space 

'Hn = where \a) and |/?) are any two elements of TCi. The space Hn is a Hilbert 

space under the inner product: 

(^1, ^2) = tr{g\ ^2) (3) 

where the trace is over a complete set of states in Tij. Mathematically, a TS, g G Tiu, can 
always be expended in terms of a basis ga/s = \o:){l3\ of Hn as 

= Y. Caf3ga/3 (4) 

A general g G Hjj may not be reducible to the "generic form" (|l]). A non-generic TS 
with the "entangled" form describes situations of a non-complete specification of the 
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conditions, that is, the final and/or initial conditions correspond to an entangled state of 
S with some other system, say S', whose degrees of freedom are traced out. In this case, 
we have two density matrix pj„ and pout, rather then two pure states as conditions. The 
conditions can be expressed as gg^t=ti = Piniti) and Q\t=t2 = Poutih)- In such circum 
stances, the occurrence of an entangled (non-generic) TS is due to the interaction of S and 
S' via the measurement device mediator, which is used to determine the conditions. Hence 
the dynamical evolution of the system is not modified (generic TS do not evolve in time 
to non-generic or vice versa). The TS of a closed system satisfies the Liouville equation: 

idtg = [H, g] (5) 

In the following we shall study the appearance of entangled two-states (^ in a dynamical 
way through the interaction of S and S'. To accommodate for this extra interaction we 
will need to modify the Liouville equation (^). 

Given by a two-state g{t) that corresponds to a pre and post selected ensemble, we can 
calculate the quantum mechanical probabilities for the result of any measurement at time 
t as follows. Let A be a Hermitian operator with a spectral expansion , A = J2 o-Pa in 
terms of projection operators Pa = \a){a\. Then, the probability to find A = a is given by 

^^^^^"'') = ^^%^^l7W^ ^6) 

Ea' \ {Pa',0{t))\^ 

Therefore, in analogy with the ordinary expression for probability, the projection of g on 
Pa, {Pa, q), can be interpreted as the TS amplitude. The absolute square of this amplitude 
is proportional to the probability. In general, this probability distribution can not be 
reduced to an expression in terms of a pure or mixed density matrix. To see this, notice 
that Equation (|) can also be written as 

-p h( (g-P"'-^'^^) _ trP^p(a) 

rob[a,t) - ^^,^^p^,^p^,^^ - j:^,tr Pa, p{a') 



where, p(a) = gPaO^ ■ Therefore, this probabihty can be expressed in terms of a density 
matrix, only when p(a) is independent of a. 

Finally, we note that if the ensemble is only pre (or post) selected, the ordinary expres- 
sion for the probability can be obtained as follows. Assuming that the final (unknown) 
measurement of some Hermitian operator K determines one of the eigenstates ipk-, the 
probability to find A = a is given by, 

Vrobi{a;t) =Y.Vrob{a)Vrob{iJk\iJin) = M^l^in)? (8) 
i.e. by the ordinary expression. In terms of the TS this yields 

trPi..(«) {i>{t),Q{t)) 
where pm = QQ^ ■ This expression is to be compared with (^. Contrary to the former case 
of a pre- and post-selection, the latter expression depends only in the initial condition. 



3 A system with a pre and post selected environment 

Consider a closed system St which is composed of the sub-systems S and Se- Let the part 
Sf, play the role of an environment £. The Hamiltonian of the total system is 

Htot = Hs + He + Hint = Hq + Hint (10) 

where Hs and He are the "free" Hamiltonians of S and respectively, and Hint is some in- 
teraction term. Given the pre- and post-selected states, \ipi{ti)) = |si)®|ei) and |V'/(^2)) = 
1^2) ^1^2), the TS in the Schrodinger representation is Qs+eit) = U {t — ti)\4'i) {ipf\W(t — t2), 
where U{t2 — ti) = exp(— i J^^ Htotdt'). Limiting out observations only to the subsystem S, 
we would like to compute the probabilities for observables of the form A = Ag ® le, where 
As operates in the Hilbert space Tis of S and Ig is a unit operator in Tig. 
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This probability can be expressed in a simple form by Eq. (|^), with q = gsj^e{t) and 
Pa = {Pa)s ® h- Obviously, since the projection operator acts only in Hs, we can trace 
over S and represent this probability in terms of a reduce TS Qs : 

Vrob{at\s2, 62, Si, ei) = I ^|^"'^^^(^^)) I (11) 



where, 



Qs = Qsit; S2, 62, Si, ei) = ^tre^s+e (12) 



The time independent normalization, = (62(^2)! exp[—ziJe(^2 — ^i)] was chosen 
for later convenience. At intermediate times S is completely described in terms of the 
reduced TS. 

Notice that at the boundaries, t = ti and t = t2, the reduced TS has a simple generic 
form: 

0s{t2) = {UUS2){S1\ = \s'){Si\ (13) 

and 

Qs{ti) = \si){s2\{lfX = \s2){s"\ (14) 

where 

ffr. (e2|?7(t2-ti)e-^-(*^-*^)|ei) 



is the 'weak value' |T^ of the evolution operator U with respect to the 'free' environment's 
pre and post-selected states. Hence, {U)w is an operator in the Hilbert space Hs- On the 
other hand, due to the interaction with the environment, at intermediate times, t G (ti, t2), 
the reduced TS is generally a non-reducible "entangled" TS: 

gs{t)=Y.Cs's"it)\s'){s"\. (16) 

This effect of "decoherence" and then "recoherence" of the reduced two-state, as expressed 
in Equations (p!3D , (jl^) , and (|16D, stands in the heart of this paper. The final post selection 
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of the environment "force's" the two-state to recohere at the final condition to a generic 
two- state. 

The effect of post-selecting the environment exists even if the sub-system S is not post- 
selected, i.e. the condition at t = ^2 is imposed only on S. In this case the probability to 
find A = a at t E (t2, ti) is given by 

Vrobia, t\r,, n, = ' ff; ^l^f (17) 

The sum above in over all possible eigenstates, {|s2)}, of an arbitrary complete set of 
operator(s) S. This probability is independent on the choice of S. 

Although, in this case, there is only one (initial) condition on S, due to the interaction 
with the pre- and post-selected environment, the sub-system S can not in general be 
described in terms of a pure or a mixed density matrix. Equation (0) can be rewritten as 

-n hf +\ \ trP„p(a) ..gs 

Vrob{a, t ra, ri, Si) = — — — -- (18) 

Ea' Pa'p{a) 

where 

p{a) = J2g{s2)PaQ\s2) (19) 

The object p corresponds to a density matrix only if it is independent of a. Intuitively, this 
happens when the condition at t = t2 on S does not "add" information. Let us examine 
this question more closely. When t — > t2, we have Qs {U)w\si){si\ = |s')(s2| and by (|19|) 



p{a,t2) = J2\s'){s2\a){a\s2){s'\ = \s'){s'\ (20) 



S2 



is independent of a. Therefore, near the final condition there is always an effective pure 
state. The initial state of the open system, |si) is mapped to a final pure state \s') by the 
"weak evolution operator" 

UM = W), (21) 
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Near the initial condition, Eq. (|TP|) yields 

p(a,ti) = {{a\{U)^{U^Ua))\s,){s,\ = C{a)\s,){si\ (22) 

The effective density matrix is proportional to a pure state, but the probability at t = ti 
depends on an unconventional normalization 

trP„p(a,ti) C{a) ^ 

Unless ilwUl, = 1, this probability may depend on the nature of the final condition on the 
environment. For example, if without the post selection we would have Vrob{a) = 1 and 
Vroh{h 7^ a) = 0, then these probability are not effected by the final post selection of £. 
But the post selection of £ does generally modify the probability in intermediat e cases as 
< Vroh{c) < 1. 

At any intermediate times, ti < t < t2, the effective density matrix (|T9|) will be a- 
dependent, and hence a complete description in terms of a unique density matrix is not 
possible. It is interesting however, that for a a limited class of observables, whose nature 
depends on the coupling with the environment, we can still construct an effective density 
matrix. To see this, let us choose the (otherwise arbitrary) set {|s2)} in Equation (0), a 
s eigenvalues of a complete set of an operators 5*^ that commute with Hint- In this case, 
for a given S2 the TS has a generic form : g{s2,t) = \s' ,t){s2,t\. Therefore, for an operator 
A = J2o,Pa which is conjugate to one of the operators Sk, we have {s2\Pa\s2) = constant. 
This implies that p{a), the effective density matrix, does not depend on a. Hence, if one 
measures only this limited class of observables, one can use the effective density matrix 
given by pdenif) = Z^sa '^^^^ density matrix is pure near the conditions at t = ti and 
t = t2, but generally corresponds to a mixed state at ^2 > ^ > ^i- 
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4 A simple example 

To exemplify these ideas we now consider a solvable model, which was used to demonstrate 
decoherence [Q, of a spin half particle (the system) coupled to spin half particles (the 
environment). Setting the free part of the Hamiltonian to zero the interaction part is taken 

as 

N 

Hi^t = Y.9kCJ.cjf^ (24) 
fc=i 

In term of the eigenstates of and (t^^\ the conditions can be expressed as 

mt = 0)) = {a\ t) + b\ i)) n(«fcl U) + Pk\ U)) = \sl)\e^) (25) 

and 

\Mt = T)) = (a'\ T) + b'\ i)) n(«fcl U) + ik)) = \s2)\e2) (26) 

The reduced TS can be derived according to Eq. (jl^), by tracing over the k = 1, ..N spins. 
The result is: 

Qs{t) = 4tt(««"x(T)| T)(T \ + bb'*x{~T)\ I 



x(o) 

+ ab'*x{2t - T)| T)(i I + ba'*x{T - 2t)\ i)(T | \ (27) 



where 

X{t') = Ui^ka'^e^^"'' + hP'^e-^'^'). (28) 



k 

At the initial and final conditions, the TS reduces to 



Qs 
and 



(t = 0) = \s,){s,\Ui{T) = -iy (a| T) + b\ i)) (a'*x(T)(T | + b'*x{-T){l |) (29) 



gs{t = T) = UUT)\s,){s,\ = ^(«X(7^)I T) + bx{-T)\ i)) ® [a'*{^ \ + b'*{i |) (30) 
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where the 'weak evolution operator' is 



x(o) 



(31) 



At intermediate times Qs{t) can not generally be reduce to a generic TS. 

Let us examine the case that only the spins (environment) are post selected. In this 
case we need to use equation ([T7|) and sum over all the final possibilities. Obviously, it is 
most convenient to sum over final eigenstates of (Jz- Hence we have two possible two-states: 

1 



and 



x(o) 



X(0) 



ax{T)\]) + hx{T-2t)\i)\®{] 



ax{2t-T)\]) + hx{-T)\[)]®{i 



(32) 



(33) 



The effective density matrix (p!9|) is in this case p(a) = Qs{})PaQ\{'\) + Qs{,[)PaQ\{,i)- Clearly, 
if we measure only ax or ay this expression reduces to pe// = |(^s(T)^KT) + ^s(i)^^(i))• 
Therefore, for these observables we have an effective density matrix: 



Peff = [ \a\\\xiT)\' + \xi2t - T)r) I T)(T 

+ l&p(lx(-T)P + |x(T-2t)p)U)ai 
+ab*[x{T)x*iT-2t) + xi2t-T)x*i-T)^\ !)(! I 
+ a*b[x{T - 2t)x*{T) + x(-T)x(2t - T)) | i)(T | 



(34) 



At the boundaries this expression reduces to 

Peff{t = T) = ^-l^(^ax(T)| ])+hx{-T)\ i)) ® {ax\T){] \+hx*{-T){[ |) (35) 

and 



Pe//(t = 0) 



2|X 



|X(T)P + |X(-T)P 



«l + i))®(a*(T |+6*a l) (36) 
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The initial and final effective density matrix corresponds to a pure state. However notice 
that the norm of the initial and final pure state is not the same. This reflects the non- 
unitarity of the 'weak evolution operator'. Hence, for a limited set of observables, we 
obtained a description in terms of a density matrix which initially decoheres and finally 
recoheres back to a pure state. It is now amusing to note that by fixing the initial and 
final states of the N spins to satisfy: |Q;fcQ;^*p = \(3k/3'j^\'^ = 1/2, (e.g. pre and post selection 
of cr^^^ = 1, {k = 1,..,N)), we can arrange that near the in itial condition, the state 
of the system is described for any observable by a pure state. In this case the system in 
intermediate time is (effectively), for some observables, in a mixed state, while for other 
observables, even a mixed state not exits. The system always 'recoheres' back to a pure 
state. 



5 Reduced two-state dynamics 

The two-state of a closed system satisfies a Liouville Equation. By focusing on a subsystem, 
and tracing over the environment's degrees of freedom we will also modify the equation 
of motion of the the reduced two-state. Some additional terms are now necessary to 
accommodate for the effect of the 'external' e nvironment. This problem is reminiscent to 
the well studied issue of environment induced decoherence. There is however a significant 
difference between the two problems. As we have seen, when the conditions correspond to 
pure states, the exact solution for the TS must be of generic (direct product) form, both, 
initially at t = ti and finall y at t = t2. Therefore, the resulting dynamical equation must 
have the non-trivial property that given any two conditions for S, it evolves an initially 
generic TS to an "entangled TS" at intermediate times, and back to a generic TS at at 
the final condition. Such a 'fine tuning' requires cushion when approximations are used 
to derive the corrections to the Liouville Equation. For example, in deriving the equation 
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of motion to the reduced density matrix, it is usually assumed that one can use the 'non- 
reversible' approximation that the density matrix can be factorized to a product of two 
density matrix of form: pdensUy = Ps ^ Pe- This simplifies considerably the computations. 
However, in our case such an approximation is invalidated since the TS can not be factorized 
in such a way at any time. In fact a naive usage of such a factorization leads to an equation 
of motion with no solutions for the two boundary condition problem. 

In the following we shall derive perturbatively the modified Liouville Equation. There- 
fore we expect our solution to be valid only in the weak coupling regime AT < 1, where A 
is the coupling constant ( Hint = Ai?/), and T = t2 — ti. For simplicity we shall assume a 
time independent Hamiltonian and that Hint is an analytic function. In the following, it 
will be most convenient to use the interaction representation. Setting ti = and t2 = T 
we define the TS in the interaction representation as 



Qint{t) 



JHot 



e{t)e 



-iHot 



(37) 



The equation of motion of the closed system is 



dtQ 



int 



[Hint] I, 



int 



(38) 



where [0]i = e*^o*[0]e-*^o*. 

Now define = |V'i(0))(V'2(^)| exp{—iHQT) = q^q (g) ^eO; which is the free {Hint = 0) 
two-state at t = 0. In terms of we have 



Qint{t) 



^-iHt^+iHot 



^-iHo{t-T)^+iH{t-T) 



(39) 



For simplicity let us assume that [Hq, Hint] = 0, hence 



(40) 



Although the exact solution gint can not be factorized, we can use (^0]) to expend it in 
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powers of = Qso x ^eo- Putting from now on ^ = Qint and XHi = [Hint]i, we have: 



g(t) = Qo- iXtHjQo - iX{T - t)goHj + 0(A^ 



(41) 



The free TS is factorizable, and we can now trace over S. Therefore, 



treg 

tre^eO 



(42) 



where {...)w stands for the weak value with respect to free environment's two-state, and is 
defined by = trO^eo/tr^eo- The last expression can be also inverted to 



Qso = Qs{t) + iXt{Hi)^Q,{t) + iX{T - t)Q,{t){Hi)^ + 0{X^) 



(43) 



Substituting ( ^T]) into the Liouville equation and tracing over £ yields 



dtQs{i) = -iX[{Hi)^, Qso] - X^ 



(Hi,tHiQso + {T -t)QsoHi 



+ 0(A=^) (44) 



Finally, we can use (|^) to reexpress the last equation in terms of git). We get 



dtQs{t) = ~iX[{Hi)^,gs{t)] 



(if7,tif7^,(t) + (T-t)^,(t)if,) J + [(if,)^,(tif,^,(t) + (T-t)^,(t)if, 



(45) 



-O(A^) (46) 



Let us consider some examples. For a generic interaction of the form: 



Hj — XQiLi 



(47) 



where the Qj's are some system variables and Li reservoir variables, we get in the free case 
{H, = H, = 0): 

dtUt) = -^X{LX[Q^, Qs] - A'A,,[g„ tQ^g^ + (T - t)QsQj\ + 0{X^) (48) 



where 



{^LlLj^yj (Z/j)^(Z/j)^ 
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(49) 



Typically, the first order is the Liouville equation with a "weak" Hamiltonian" . The 
second order corrections, are proportional to the "weak uncertainty" Aj^. Higher order 
may be easily computed, but become very cumbersome. It is straightforward to rewrite 
(^) to the case that Li and Qi are not constants of motion, or to any other polynomial 
interaction. 

Simplifying the interaction even further, we set Qi = a^, Li = = L and = Qi = 
for i ^ 1. This corresponds to a spin half subsystem which interacts with the z component 
of the angular momentum of the environment. Equation (^) reduces to 

dtQs = -tXL^a,, Qs] - 2A2AL^(2t - T)(^, - a.^.a,) + OiX^) (50) 

where AL^ = (L^)^ - 

We can easily verify that for every two initial and final conditions for S, there exists 
an appropriate solution. It is only the second order term that can induce transition from 
generic to non-generic (entangled) two-state. In terms of the notation g-^^ = \ |)(| |, etc, 
the general solution of Eq. (|50|) is 

= ^TTO, hd^) = hio (51) 



Quit) = exp 



QUO (52) 



-i2XL^t - AX^AL^if - Tt) 

Q^^it) = exp +i2XL^t - AX'^AL^if - Tt) Qi^q (53) 

Clearly, due to the factor t"^ — Tt, the second order contributions vanishes on at the con- 
ditions. By substituting XL = J^Qk^z, it can be verified that this solution agrees up to 
corrections of order 0{X'^) with the exact solution given by equation (P?]). 

Due to the continues interaction with each of the spins in the latter problem, the validity 
of Equation ([50| ) is limited by the constraint T < 1/X. We shall now compare this system 
to the other extreme case, in which the subsystem interacts with each of the particles of 
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the environment separately, and only for a very short time At = r, such that rA << 1. 
In this way the weak coupling condition is satisfied, and our modified Liouville Eq. can 
be applied also for long tim es. Let the environment be composed of non-interacting 
particles. The interaction Hamiltonian for this case is given by 

N 

Hr = \J2U{t)H^ (54) 

n=l 

where fn{t) = 6{t ~ nr) — 9{t — {n + 1)t) with 9{t) as the step function is nonzero only 
for t G (nr, (n + 1)t). Hn is the interaction of S with the nth particle. Let us further 
assume that can be regarded as (or i s) constant during the interaction times r. For 
Hn = cr^Lnz = crLn we get 

dtQs = -i\^fn{t)Lnw[<J, Qs] - A^^/„A„„(2t - (2n + 1)t){qs - aQsCr) 

n n 

m=n—l m=N 

- fn^nm{nT)[a, aQs] - ^ /„A„„(A^ - n - l)r[cr, gs(r] (55) 

n,m=l n,m=n+l 

where A„m = {LnLm)w ~ {Ln)w{Ln)w If the initial and final states of the environment are 
given by a product state, Hfc ®|efc) of the particles, there are no correlations between the 
weak values different particles in the reservoir and Anm = ((-^n)^ ~ {Lnwy)Snm- Therefore, 
in this case the two last terms on the right hand side of equation (^) vanish. Integrating 
(p5|) we see that after each "step", when the interaction with the n'th particle in the 
environment is completed, the accumulated contribution of the second term drops to zero. 
The TS remains 'pure' up to fluctuations of order O(A^r^). In this sense, we can say that 
the post-selection of the environment prevents decoherence of the subsystem, 
t a time scale T. 
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